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Abstract We considered the effects of convection on 
the radiatively inefficient accretion fiows (RIAF) in the 
presence of resistivity and toroidai magnetic field. We 
discussed the effects of convection on transports of an- 
gular momentum and energy. We established two cases 
for the resistive and magnetized RIAFs with convec- 
tion: assuming the convection parameter as a free pa- 
rameter and using mixing-length theory to calculate 
convection parameter. A self-similar method was used 
to solve the integrated equations that govern the be- 
havior of the presented model. The solutions showed 
that the accretion and rotational velocities decrease by 
adding the convection parameter, while the sound speed 
increases. Moreover, by using mixing-length theory to 
calculate convection parameter, we found that the con- 
vection can be important in RIAFs with magnetic field 
and resistivity. 

Keywords accretion, accretion discs, convection, 
magnetohydrodynamics: MHD 



1 Introduction 

The existence of radiatively inefficient accretion flows 
(RIAFs) have been confirmed in low-luminosity state 
of X-ray binaries and nuclei of galaxies (Narayan et al. 
1996; Esin et al. 1997; Di Matteo et al. 2003; Yuan 
et al. 2003). It was understood that RIAFs are likely 
to be convectively unstable in the radial direction due 
to the inward increase of the entropy of accreting gas 
(Narayan & Yi 1994). Moreover, hydrodynamical and 
magctohydrodynamical simulations of low- viscosity RI- 
AFs have confirmed these flows are convectively unsta- 
ble ( e. g. Igumenshchev et al. 1996; Stone et al. 1999; 
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Machida et al. 2001; Hawley & Balbus 2002; McKin- 
ncy & Gammic 2002; Igumenshchev et al. 2003). Self- 
similar or global solutions for convection-dominated ac- 
cretion flows (CDAFs) were presented by several au- 
thors (e. g. Narayan et al. 2000; Quataert & Gruzinov 
2000; Abramowicz et al. 2002; Lu et al. 2004; Zhang 
& Dai 2008). 

Igumenshchev et al. (2003) studied the resistive 
MHD simulations of RIAFs onto black holes. They as- 
sumed two cases for the geometry of the injected mag- 
netic field: pure toroidal field and pure poloidal field. 
They found that in the case of pure toroidal magnetic 
field, the accreting gas forms a nearly axisymmctric, ge- 
ometrically thick, turbulent accretion disc. Moreover, 
their solutions represented that the flow resembles in 
many respects CDAFs found in previous numerical and 
analytical investigations of viscous hydrodynamic flows. 
Zhang & Dai (2008) investigated the effect of magnetic 
field on RIAFs with convection by a scmi-analytically 
method. By exploit of a-prescription for viscosity and 
convection, they used two methods to study of mag- 
netized flows with convection, i.e. they take the con- 
vective coefficient a c as a free parameter to discuss the 
effects of convection for simplicity. They also estab- 
lished the a c -a relation for magnetized flows using the 
mixing-length theory and compare this relation with 
the non-magnetized case. They found that the mag- 
netic field makes the a c -a relation be distinct from that 
of non-magnetized flows. 

Since the importance of toroidal magnetic field and 
resistivity in accretion flows have been confirmed ob- 
servationally (see Faghei 2011 and references therein), 
Faghei (2011) considered the steady, radially self- 
similar solutions of accretion flows in the presence of 
the toroidal magnetic field and the resistivity. How- 
ever, he ignored the effects of convection in his model. 
Generally semi-analytical studies of magnetized CDAFs 
are related to non-resistive magnetized CDAFs (e. g. 
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Zhang & Dai 2008) and the resistive and magnetized 
CDAF was studied in MHD simulations (e. g. Igu- 
menshchev & Narayan 2002; Hawley & Balbus 2002; 
Igumcnshchev et al. 2003). Thus, it will be interest- 
ing to study the effects of resistivity on RIAFs with 
convection. Here, we adopt the presented solutions 
by Narayan ct al. (2000) and Faghei (2011). Similar to 
Narayan et al. (2000), we will discuss the effects of con- 
vection on angular momentum and energy equations. 
The paper is organized as follow. In section 2, the ba- 
sic equations of constructing a model for quasi-spherical 
magnetized RIAFs with convection will be defined. In 
section 3, a self-similar method for solving equations 
which govern the behavior of the accreting gas was uti- 
lized. The summary of the model will appear in section 
4. 



2 Basic Equations 

Analytical theory of CDAF is based on a self-similar 
solution of a simplified set of equations describing RI- 
AFs. We adopted the presented solutions by Narayan 
et al. (2000) and Faghei (2011). By using spherical co- 
ordinate (r, 9, ip) centered on a accreting object, let 
us consider stationary, axisymmctric, quasi-spherical 
equations describing an accretion flow onto the black 
hole of mass M. For the sake of simplicity, the general- 
rclativistic effect has been neglected and the gravita- 
tional force on a fluid is characterized by Newtonian 
potential of a point mass, ip = —GM/r. As magnetic 
fields, we consider only toroidal fields, B v . 

Under these assumptions, the continuity equation 
with mass loss is 
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where p, v r and p are the density, the accretion ve- 
locity (v r < 0), and the mass- loss per unit volume, 
respectively. 

The radial momentum equation is 
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where c s is sound speed, which is defined as c 2 = p gas / p, 
with being p gas as the gas pressure, ft is the angular 
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velocity, Qk[= (GM/r 3 ) ] is the Keplerian angular 
velocity, and ca is the alfven speed, which is defined 
as c 2 A = B^/Airp = 2p mag /p, with being p mag as the 
magnetic pressure. The ram-pressure term v r dv r /dr 
and last two terms due to the magnetic field in this 



equation were ignored in the self-similar CDAF model 
of Narayan et al. (2000), while we include them here in 
order to consider their effects. 

The angular momentum equations can be written 
in the form of the balance of advection and diffusion 
transport terms (Narayan et al. 2000), 
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where the two terms of right hand side represent the 
angular momentum transport by viscosity and convec- 
tion. Here, v is the kinematic viscosity coefficient, v c is 
the convective diffusion coefficient, and g is the param- 
eter to determine the condition of convective angular 
momentum transport. When g = 1, the flux of angular 
momentum due to convection is 
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The above equation implies that the convective angu- 
lar momentum flux is oriented down the angular ve- 
locity gradient. For a quasi-Kcplcrian angular velocity, 
Q oc r~ 3 / 2 , angular momentum is transported outward. 
When g = —1/3, the convective angular momentum 
flux can be written as 



Jc = -v c pr 



d (Qr 2 ) 
dr 



(5) 



This equation represents that the convective angular 
momentum flux is oriented down the specific angular 
momentum gradient. For a quasi-Keplerian angular ve- 
locity, fl oc r~ 3 / 2 , angular momentum is transported 
inward. Generally, convection transports angular mo- 
mentum inward (or outward) for g < (or > 0), and 
the specific case g = corresponds to zero angular mo- 
mentum transport (Narayan et al. 2000). 

In this paper, we assume the kinematic coefficient of 
viscosity and the magnetic diffusivity due to turbulence 
in the accretion flow. So, we use these parameters in 
analogy to the a-prcscription of Shakura & Sunyaev 
(1973) for the turbulent, 
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where P m is the magnetic Prandtl number of the turbu- 
lence, which assumed to be a constant less than unity, 
rj is the magnetic diffusivity, and a is a free parameter 
less than unity. For the convective diffusion coefficient, 
v Cl we adopt the assumptions of Narayan et al. (2000) 
and Lu et al. (2004) that all transport phenomena due 
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to convection have the same diffusion coefficient, which 
is defined as 



T 2 

M 
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where Lm is the characteristic mixing length and N e f f 
is the effective frequency of convective blobs. The char- 
acteristic mixing length Lm in terms of the pressure 
scale height, H p , can be written as 



Lm — 2 x ^ImH p , H„ — — 



dr 



d\np gas 



(8) 



where Im is the dimensionless mixing-length parameter 
and its amount is estimated to be equal to y/2 in ADAFs 
(Narayan et al. 2000; Lu et al. 2004). the effective 
frequency of convective blobs, N e ff, is given by 



N 2 ff = N 2 + K 2 , 



(9) 



where N is Brunt- Vaisala frequency, which is defined 
as 
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and k is epicyclic frequency, which is defined as 
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For a non-Keplerian flows k 7^ 57, while for a quasi- 
Keplerian (0 oc r~ 3 / 2 ), n = f2 (Narayan et al. 2000; 
Lu et al. 2004). Convection appears in flows with 
N 2 fj < 0. We also write the convective diffusion coef- 
ficient in the form similar to usual viscosity of Shakura 
& Sunyaev (1973), 



-V • F conv , with being F conv [= -pv c Tds/dr] as the 
outward energy flux due to convection. For the right 
hand side of the energy equation, we can write 
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where Q a dv is the advective transport of energy, and 
/[= 1 — Qrad/Qdiss] is the advection parameter. The 
parameter / measures the degree to which the flow is 
advection-dominated (Narayan & Yi 1994). The dissi- 
pative heating rate can be written as 
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where the right-hand side terms are heating rate due 
to viscosity, convection, and resistivity, respectively. In 
above equation, J [= V x B] is the current density, with 
being B as the magnetic field. 

Finally, the induction equation with creation/escape 
of magnetic field can be written as 
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where B v is the toroidal component of magnetic field 
and B v is the field escaping/creating rate due to a mag- 
netic instability or dynamo effect. This induction equa- 
tion is rewritten as 
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where x is the ratio of the magnetic pressure to the gas 
pressure, which is defined by 
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where a c is a dimensionless coefficient that describes 
the strength of convective diffusion. The a c coefficient 
can be obtained by equations (8) and (13) 
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The energy equation is 
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where T is the temperature, s is the specific entropy, 7 
is the ratio of specific heats, Qdiss is dissipative heating 
rate, Q ra d is the radiative cooling rate, and Q con v = 
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3 Self-Similar Solutions 

We seek self-similar solutions in the following form (e.g. 
Narayan & Yi 1994; Akizuki & Fukue 2006) 
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where ci, C2, and c 3 are dimcnsionlcss constant to be 
determined. We use a power-law relation for density 



p{r) = par* 



(24) 



where po and A are constant. Using equations (20)- 
(24), the mass-loss rate and the magnetic field escap- 
ing/creating rate can be written as 



p(r)=p r^/\ 
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Above equations express for A = —3/2, there is no mass 
loss, while for A > —3/2 mass loss (wind) exists. 



4 Results 

Here, similar to Zhang & Dai (2008), we will study 
the presence of convection in two cases: a c as a free 
parameter and a c as a variable. 



B v (r) = B r' 



(26) 



where po and Bq are constant. Since we have not ap- 
plied the effects of wind in the momentum and energy 
equations, we will assume a no wind case, p = and 



_11 / 4 , which implies 



A = —3/2. In this case, B v oc r 
that creation/escape of magnetic field increases with 
approaching to central object. This property is qual- 
itatively consistent with previous studies of accretion 
flows (Machida et al. 2006; Oda et al. 2007; Faghei & 
Mollatayefeh 2012). 

Using the self-similar solutions in the continuity, ra- 
dial momentum, angular momentum, convection pa- 
rameter, energy, and induction equations [(l)-(3), (13), 
(14), and (18)], we can obtain the following relations: 
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4.1 Case 1: a c as a free parameter 

In this case, we take the convective coefficient a c as a 
free parameter to discuss the effects of convection for 
simplicity. Examples of such solutions are presented in 
Figures 1 and 2. 

In Figure 1 , the self-similar coefficients c\ , ci , and c 3 
are shown as functions of the parameter \. By adding 
the parameter \ which indicates the role of magnetic 
filed on the dynamics of accretion discs, we see the co- 
efficients of radial and rotational velocities and sound 
speed decrease. This properties are qualitatively con- 
sistent with results of Faghei (2011). In Figure 1, we 
also studied the effect of convection parameter a c on 
the physical variables. The value of a c measures the 
strength of convective viscosity and a larger a c denotes 
a stronger turbulence due to convection. Figure 1 im- 
plies that for non-zero a c , the radial infall velocity is 
lower than the standard ADAF solution and for larger 
ct c this reduction of radial infall velocity is more evi- 
dent. It can be due to decrease of efficiency of angular 
momentum transport by adding the convection param- 
eter a c (sec equation 29). The profiles of angular ve- 
locity show that it decreases with the magnitude of a C: 
while the sound speed increases. These properties are 
in accord with results of Zhang & Dai (2008). 

In Figure 2, the physical variables are shown as func- 
tions of parameter \ f° r several values of magnetic 
Prandtl number. Since inverse of magnetic Prandtl 
number is proportional to magnetic diffusivity, P m oc 
rj^ 1 . Thus, reduce of magnetic Prandtl number denotes 
to increase of resistivity of the fluid. The solutions in 
Figure 2 imply that the accretion velocity and the sound 
speed both increase with the magnitude of resistivity, 
while the rotational velocity decreases. These proper- 
ties qualitatively confirm the results of Faghei (2011). 





4.2 Case 2: a c as a variable 

Here, we calculate the dimensionlcss coefficient a c by 
using the mixing-length theory. Because we used a 
steady self-similar method to derive a c , it becomes a 
constant throughout of the accreting gas. However, it 
is a function of position and time (e. g. Lu et al. 2004). 
The amount of convection parameter a c is calculated by 
equation (31). Using this equation and equation (28)- 
(30), we can obtain the behavior of physical quantities 
in the presence of convection. Such solutions are shown 
in Figure 3. 

In Figure 3, the coefficients c\, C2, C3, and convec- 
tion parameter a c are shown as functions of the degree 
of magnetic pressure. Similar to case 1, the accretion 
and rotational velocities, and sound speed decrease by 



adding the parameter \- While, the convection param- 
eter a c increases for stronger toroidal magnetic field. 
This property is qualitatively consistent with result of 
Zhang & Dai (2008). In Figure 3, the physical variables 
are also studied for several values of magnetic Prandtl 
number. The profiles of convection parameter a c imply 
that it increases by adding the magnetic diffusivity. As 
for non-zero magnetic diffusivity, a c is larger than the 
standard CDAF solution and for larger magnetic diffu- 
sivity this increase of convection parameter a c is more 
evident. 

5 Summary and Discussion 

The observational features of low-luminosity state of X- 
ray binaries and nuclei of galaxies can be successfully 




Fig. 3 Physical variables as functions of \ f° r several values of magnetic Prandtl number. The input parameters are set 
to a = 0.5, 7 = 1.5, f = 1, I = v2j <? = — 1/3, and A = —3/2. The solid, dashed, and dotted lines represent P m = oo, 1.0, 
and 0.5, respectively. 



explained by the models of radiatively inefficient ac- 
cretion flow (RIAF). The importance of convection in 
RIAFs was realized by semi-analytical and direct nu- 
merical simulation (e. g. Narayan et al. 2000; Igumen- 
shchev et al. 2003). 

In this research, we considered the effects of convec- 
tion on the presented model of Faghei (2011). Simi- 
lar to Narayan et al. (2000), we assumed the convec- 
tion affects on transports of angular momentum and 
energy. Using a radially self-similar approach, we stud- 
ied the effects of convection on the model for several 
values of magnetic field and resistivity. The solutions 
showed that the accretion and rotational velocities, and 
sound speed decrease for stronger magnetic filed. More- 
over, we found that the accretion velocity and sound 
speed increase with the magnitude of the resistivity, 
while the rotational velocity decreased. These proper- 
ties are qualitatively consistent with results of Faghei 
(2011). We studied the effects of convection on a re- 
sistive and magnetized RIAF in two cases: assuming 
the convection parameter MS (X free parameter and us- 
ing mixing length theory to calculate the convection 
parameter. In the first case, we found that by adding 
the convection parameter, the radial and rotational ve- 
locities decrease and the sound speed increases. In the 



second case, we found that the convection parameter 
increases by adding the magnetic filed and resistivity. 
These properties are in many aspects in accord with 
results of Zhang & Dai (2008). 

The present model have some limitations that can 
be modified in the future works. For example, the lati- 
tudinal dependence of physical variables have been ig- 
nored in this paper. While, two-dimensional and three- 
dimensional MHD simulations of RIAFs show that the 
disc geometry strongly depends on magnetic field con- 
figuration (e. g. Igumcnshchev et al. 2003). Thus, the 
study of present model in two/three dimensions can be 
an interesting subject for future research. Moreover, 
it has been understood the magnetic field can change 
the criterion for convective instability (e. g. Balbus & 
Hawley 2002). While, we igonred the effects of mag- 
netic field on the instability criterion. Thus, the pre- 
sented criterion in this paper can be modified in the 
future research. 
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